Novel reconstruction technique for New Physics processes with initial state radiation 
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At hadron colliders, the production of heavy new particles is associated with additional quarks and 
gluons with significant transeverse momentum. The additional jets complicates the reconstruction 
of new particle masses. Taking gluino pair production and decay at the Large Hadron Collider as an 
example, we develop a novel technique to reduce these effects, and to reconstruct a clear kinematical 
endpoint for the gluino decay products. 



Although the Standard Model very well describes the 
interactions of elementary particles, the Higgs boson, the 
particle responsible for the electroweak symmetry break- 
ing, is associated with a hierarchy problem. Many mod- 
els have been proposed to solve this problem, some of 
which can also account for the dark matter inferred by 
astronomical observations. Examples of such models are 
Supersymmetric models with R parity and Little Higgs 
models with T parity. These models predict the existence 
of new colored particles around 1 TeV, as well as a stable 
lightest new particle. 

In supersymmetric models, squarks and gluinos will be 
copiously produced in the ATLAS and CMS experiments 
at the LHC By looking at kinematical distributions 
of the decay products, the masses of squarks and gluinos 
can be reconstructed. However, since the production and 
decay processes involve particles charged under QCD, ini- 
tial and final state QCD radiation will complicate this 
reconstruction 0. 

The problem of initial state radiation (ISR) is generic 
for all new physics signatures involving multiple hard jets 
in the final state. It reduces the precision on mass de- 
terminations in hadronic channels. Understanding of jet 
emission associated with hard processes is becoming one 
of the central issues of the LHC phenomenology. 

In Supersymmetric models, the effect of initial state 
radiation is especially severe for the pair production of 
gluinos. The lowest order process pp — > gg is dominated 
by gg — * gg- When both gluinos decay into qqxl, it 
is possible to determine nig and m^o from the jet dis- 
tributions. In particular, it was recently realized that 
the gluino and LSP masses can be simultaneously recon- 
structed by looking at the end point of Mti distributions 
[H, H[, if we can correctly identify two jet pairs that arise 
from gluino decays. However, by picking up jets com- 
ing from initial state radiation, the end point is signifi- 
cantly smeared, leading to increased uncertainty in the 
end point determination. 

The aim of this paper is to propose a new jet selec- 
tion method to reduce the misreconstruction due to the 
initial state radiation. Systematic studies of the effect of 
additional QCD radiation has recently become possible 
thanks to Monte Carlo (MC) simulations which include 



the effect of hard parton emission frorn inital state quarks 
and gluons in new physics processes [f|. 

Processes involving initial state radiation can be simu- 
lated using parton shower (PS) MCs such as PYTHIA 
or HERWIG. They generate soft and collinear parton 
shower emissions in association with hard two-to-two pro- 
cesses. Hard parton emissions, on the other hand, are 
not correctly described by this approximation. There, 
matix element (ME) calculations are required to cor- 
rectly predict distributions. Various matching techniques 
have been developed to remove double counting between 
PS emission and ME contributions and account for Su- 
dakov suppression effects @. These techniques are im- 
portant to compare the differential distributions of the 
hard jets with theoretical predictions. Detailed com- 
parisons of various event generators with PS-ME match- 
ing have been performed [7j, showing that the matching 
schemes used in Sherpa, MadGraph, ALPGEN, and ARI- 
ADNE are consistent each other, although the algorithms 
are slightly different. 

For SUSY processes, PS and ME matching is avail- 
able in MadGraph/MadEvent Q matched with Pythia 
We here use it to study the effect of initial state 
radiation in sparticle mass reconstruction at the LHC. 
We have generated sparticle production events with up 
to one additional quark/gluon. The matching scheme 
used is called "A^-jet MLM" scheme, described in de- 
tail in [5(. In this scheme, the final state partons in an 
event are clusted according to the kx jet algorithm to find 
the equivalent parton shower history of the event. The 
smallest Ut value is restricted to be above a cut off scale 
Q^f. After showering, the final state partons are clus- 
tered into jets using the fcx jet algorithm with a cutoff 
scale Q m atch > Qcut • The event is rejected unless each jet 
is matched to a parton, except for the highest multiplic- 
ity sample. Events from the pp — » XX' contribution to 
the matched sample, where X and X' is either a gluino 
or a squark, are here referred to as 'exclusive' events, 
while events from the highest multi-jet sample (in our 
case XXj), which contains hard emissions as well as re- 
solvable PS effects, are called 'inclusive'. In this paper we 
take Q match =60 GeV and Q^f = 40 GeV, so 'inclusive' 
events have additional hard radiation with pt >60 GeV. 
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For demonstration purposes, we focus in this paper 
on the study of gluino pair production. For simplicity 
we force the gluinos to decay into uuxi- In this con- 
text, the exclusive sample is events with four partons 
+ 2 LSP at parton level, and the inclusive sample is 
events with five partons + 2 LSP. Squark masses are 
taken to be high enough so that their production can be 
ignored. We take an MSSM point with m g = 685 GeV 
and mq = 1426 GeV, and m x o = 101.7 GeV. We gen- 
erate 1.21 x 10 5 events. After removing squark gluino 
associated production (where an on shell squark decays 
into gq to generate the ggq final state), 1.07 x 10 5 events 
remains. The cross section is 2.5 pb, so the number of 
generated events corresponds to 40 fb _1 . 

The generated events are then simulated by the toy de- 
tector simulator AcerDET [lOj with the jet reconstruc- 
tion tool Fastjet [ll|. In AcerDET, the phase space is 
divided into cells with (A?y, A(j>) — (0.1, 0.1), and the mo- 
menta of hadrons, electrons, and photons passing within 
one cell are summed to imitate the energy deposit in a 
calorimeter cell. The cell energy deposits are interfaced 
to Fastjet for jet reconstruction as massless particle mo- 
menta. The momenta of the reconstructed jets are then 
smeared as p sme ar = (l+^)Pjct, with the energy resolution 
S = 0.5(l)/-y/-Ej c t assumed in the barrel (forward) direc- 
tion with |fy| < (>)3.2. We use the Cambridge- Achen 
algorithm with R = 0.4 in this paper. 1 The missing 
Pt was calculated from the reconstructed objects after 
smearing. 

In Fig. [T]a), we show the Mt2 distribution calculated 
using the four hardest jets. The Mti observable is cal- 
culated from two visible momenta p\ ls , p\ ls , a test LSP 
mass m x st , and test LSP momenta p\ x andp2x satisfying 
the constraint pi x + p2 X = p™ lss as follows, 

Mt2 = min 

[max (MT(p v Apf x \ mf*), M T {p? tI % x , mf*)] . 

(1) 

If both p\ ls and are the momenta of the sum of a 
visible gluino decay products and m^ ost is taken as the 
LSP mass, the Mt2 endpoint should coincide with trig. 
It has also recently been pointed out that the Mt2 end 
point, taken as a function of the test mass m x st , shows a 
kink at the true LSP mass. Therefore, the LSP mass and 
gluino mass can be determined simultaniously. These 
ideas have got significant attention in the literature, and 
various extensions are beeing studied. 

Experimentally, we cannot know from which parent 
particle a jet arises. In Fig. [T]a), we have defined p v - ls as 
follows: 



1 In principle, the jet resolution should depend on the reconstruc- 
tion algorithm and R. In addition, energy smearing and momen- 
tum smearing may not be the same. The smearing we introduced 
in this work is only for illustrative purposes. 




UT2(GeV) UT2(GeV) for the highest four pt jets. r>50=4 



FIG. 1: a) (left) The rriT2 distribution calculated from the 
four highest pr jets, using all events, b) (right) The Mt2 
distribution of events with exactly 4 jets with pr > 50 GeV. 

1. We first take the two highest px jet momenta p\ 
and P2 as seeds. 

2. We then calculate Mt2 for the combinations 1) 
(pf s ,P V 2 is ) = (Pi + Pa, P2 + Pi) and 2) (p! + 

1 ( 2) 

P4i P2 + P3) giving M T2 and take the minimum, 
M T2 = min(M^ 2 , M T2 ). 

In the Fig. [TJ a) we show the Mt2 distribution of the 
four highest px jets paired as above for the events with 
"so = (number of jets with pt > 50 GeV) > 4 as the 
solid line. The distribution does not show a clear end 
point. This is due to the effects of the initial state ra- 
diation. To show this, we plot the distributions of the 
exclusive and inclusive samples separately in the same 
figure. A dashed line shows the distribution of exclusive 
events, i.e., events in the matched sample which have 
no resolvable QCD radiation above 60 GeV (4-parton 
events). The dotted line shows the distribution for the 
inclusive sample, i.e., the contribution from five parton 
events. The exclusive sample has an end point at the cor- 
rect gluino mass, since it does not contain hard additional 
jet activity besides that coming from gluino decays. The 
ratio iV(5 parton) /N(A parton) is here 1.4. The fraction 
of the events with additional partons after the matching 
is larger for gluino pair production compared with squark 
pair production. For example, at the SPS la point, where 
the gluino mass is 595 GeV and squark mass is around 
530 GeV, The ratio N(XXj)/N(XX) is 1.47 for X = g 
and 0.81 for X = Q [= u, d, c, s and its charge conju- 
gates) for a 60 GeV jet resolution scale. 

There are several reasons why the smearing of the end 
point due to initial state radiation is here relatively large. 
In the simulation, the gluino is forced to decay into a 
three body final state. The typical pt of the partons 
from gluino decay is therefore ~ rrig/3. On the other 
hand, the pr of the initial state radiation is not small, 
in average O(100) GeV, because the produced gluino is 
heavy. The pt distribution is shown in Fig. [2] a). The 
Pt of the additional parton is large in average, and often 
larger than one of the partons from gluino decay, as can 
be seen in Fig. [2] b) where the position of of the initial 
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FIG. 2: a) (left pr distribution of the additinal parton for 
PP 993- The dotted(dashed) line shows the distribution for 
j = g(q). b) (right) pr order of the ISR parton among the 
five parton of the inclusive sample. 
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FIG. 3: a) (left) Distribution of M™ 2 in at jet level for the 
events with 7150 > 5 and i m i n > 3. b) (right) parton level 
M™ distribution for the 5 parton sample. 

state parton among all partons, ordered in pt, is plotted. 
The probablity that ISR parton is the 5th, softest, parton 
is only about 22% of all five parton events. 

One may recover the clean end point by requiring ex- 
actly 4 jets with pr > 50 GeV in the final state (Fig. [Tj 
b). However, this selection is not practical for general 
MSSM model points. The reason is that we expect the 
decay branching ratio into heavier neutralino or chargino, 
9 -* xU^iX + qq\ to be large, where Xi,xt further de- 
cays into jets and lcptons. The branching ratio for both 
of the gluinos to decay into 2 jets and LSP may be small, 
in which case this cut would reduce the statistics signifi- 
cantly. 

A better solution is obtained by taking into account the 
existence of additional ISR jets in the analysis. Given the 
high probability to have ISR jets, we should regard the 
process we are interested in as a five-jet system rather 
than the four-jet system expected at the lowest order. 
We therefore propose the reconstruction of a five jet dis- 
tribution, rather than four jet distribution hitherto con- 
sidered. 

For this purpose, we define Mrzii) [i = l...,5) ; where 
Mt2{i) is calculated from the five highest pt jets, exclud- 
ing the i— th highest pr jet. 

M T2 (i) = M T2 (pi, ...,Pi~i,p i+ i, ...p 5 ) (2) 

We find that the problem arising from the ISR is sig- 
nificantly reduced if we look at the M™ 2 ln distribution, 
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FIG. 4: The test mass dependence of the M™ 2 m end point, 
where 

M^ 2 in = min (Mraf*)). (3) 

i=l,..5 

At the parton level, M™ 2 m contains at least one correct 
parton combination, and therefore M™ 2 ln < M^ d . The 
jet level distribution is shown in Fig. [3] a) for events with 
7150 > 5 and z m i n > 3, where z m in satifsy Mx2(Wn) — 
Af™ 2 m . Note that the events in Fig. [3] a) are statistically 
independent from those in Fig.[TJb). 

In the figure, we do not include the events with i m i n = 
1, 2, because My 2 (l), and My 2 (2) tend to be softer than 
the others since the removed jet has a high px- The 
distribution of i m - nl = 1,2 appears to be smeared and 
the end point is not well determined. The parton level 
M™ 2 m distribution for the five parton sample is given in 
Fig. SI b). 

We fit a f(x) to the M™ 2 ln distributions, where f(x) = 
G(x - M cnd )[ ai (x - M ond ) + b]+0(M cnd - x)[a 2 {x - 
M end ) + b] to see if the end points are recovered cor- 
rectly. The fitted end point Af™ d is 692.3 ± 1.2 GeV at 
parton level for 5 parton events with i5,; n on > 3. This 
should be compared with M£ 2 d for the parton level Mt 2 
distribution calculated using the partons coming from 
the g decay, 694.1 ± 0.5 GeV. The M^ d end point at 
jet level is 691. 5± 3.9 GeV for events with nso > 5 and 
«min > 3. The M£ 2 d calculated for 7150 — 4 events is 692.3 
± 2.4 GeV. These values are consistent with the input 
gluino mass rrig = 685 GeV. The central value depends 
slightly on the fit region, and careful study is needed to 
identify systematical errors. As a cross check, we also 
study the distributions with fey and anti-fey jet recon- 
struction algorithms[12] and find the results are consis- 
tent with these values. We have also checked the test 
mass dependence on the end point. In Fig. [4] the bars 
show the fitted M™ 2 m end points and the error for the 
events 7150 > 5 and 7 m j n > 3. The end points are close to 
the expected end point values shown in the solid lines. 

Several comments are in order. First, to study longer 
cascade decay chains, we have to consider processes with 
more than four final state jets. For example, when g de- 
cays into heavier inos, there may be additional jets com- 
ing from the neutralino and chargino decays. With GUT 
relations for the gaugino masses, jets coming from the end 
of the cascade decay chains are softer than those coming 
directly from the gluino decay. Our assumption that one 
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FIG. 5: (a) (left) r\ distribution of ISR quark and gluon of 
ggj production (solid line) for those with pr > 100 GeV. The 
dashed line is for ISR quark distribution and dotted line is for 
ISR gluon. b) (right) Correlation between |»j| of ISR parton 
and \rj\ for the jets that gives Af™ 2 m - 

of the five highest px jets arises from ISR is therefore 
reasonable. In such a case, one may use the inclusive 
definition of Mt2 proposed in [H, [H[. In this approach, 
p v is is defined using all jets and leptons in the final state 
so that they satisfy p^J s = JYpf 5 , p^l = JYpf' where 
p^ and p' 2 ' are jets or lepton momenta which satisfy 
d(p^,p«) < dipiZp^), d(pW, p f) < d(p«,pf), with 
d being some distance measure. The inclusive Mt2 can 
also be used to determine squark and gluino masses when 
niq > nig. In the inclusive definition of Mt2, we may 
again remove one of the leading five jets, use hemisphere 
reconstruction to define p V i S , and then calculate M™ 2 m . 
This method should be useful to reduce the contamina- 
tion from squark-gluino. 

Second, contrary to naive intuition, the additonal ISR 
jet cannot be removed by excluding jets with high tj or 
low pt from the kinematical reconstructions. We have 
seen already that the average pr of the that additional 
parton is rather high. In addition, the ISR jets are cen- 



tral. In Fig. [5] a) , we show the rj distribution of the 
additinal parton for ggj. We see that gluino ISR is al- 
most central, while quark ISR is rather forward. How- 
ever, they tend to be at high energy, as can be seen in 
Fig. [2] a) . For further details we refer to . 

Finally, the proposed method makes it possible to se- 
lect ISR jets, by requiring additional cuts to Af™ 2 m . For 
an event near the Mt2 end point, the removed jet has 
a higher probability to be the ISR jet. The probability 
that a different jet combination is correct is small, be- 
cause the correct value has to be in the narrow range 
M^ 2 in < M^ c < Mf£ d . To check this, we study the 
nature of the removed parton that gives M™ 2 ln in parton 
level. Among the 5 parton events generated, only 29% 
of the partons that give the M™ 2 ln is the ISR parton, if 
no restriction is applied to M™ 2 m . This fraction increases 
to 44% for events with M™ n > 500 GeV, 29% of total 
events. In Fig. b), we show a 2-dimentional plot where 
the x-axis is the \n\ of the ISR parton and y-axis is \rj\ 
of the jet that gives M™ 2 m . The correlation is especially 
good for \r]\ > 2, roughly 65% for the forward jets that 
match correctly to the ISR parton within |A?/| < 1. This 
is because the jets from gluino decay mostly goes to the 
central regions. This shows it is possible to study forward 
ISR jet distributions associated with the hard process. 

In this paper we have seen that ISR is an important 
feature in g production at the LHC, and we have devel- 
oped a method to reduce the effect of ISR production on 
the gluino mass determination. This method can also be 
used to identify initial state radiation jets. The method 
can be applied for any new physics processes. The ap- 
plication of this method to other SUSY processes as well 
as other models for new physics, and to the correspond- 
ing Standard Model backgrounds, will be discussed in 
following publications [l5[ ■ 



[1] G. Aad et al [The ATLAS Collaboration], arXiv: 
0901.0512 [Unknown]. G. L. Bayatian et al [CMS Col- 
laboration], J. Phys. G 34, 995 (2007). 

[2] T. Plehn, D. Rainwater and P. Skands, Phys. Lett. B 
64 5, 217 (2007 ). A. Papaefstathiou and B. Webber, 
larXiv:0903.20T3l [hep-ph] . 

[3] A. Barr, C. Lester and P. Stephens, J. Phys. G 29 (2003) 
2343. 

[4] W. S. Cho, K. Choi, Y. G. Kim and C. B. Park, Phys. 

Rev. Lett. 100 (2008) 171801. 
[5] J. Alwall, S. de Visscher and F. Maltoni, JHEP 0902 

(2009) 017. 

[6] S. Catani, F. Krauss, R. Kuhn and B. R. Webber, JHEP 
0111 (2001) 063. M. L. Mangano, M. Moretti, F. Pic- 
cinini, R. Pittau and A. D. Polosa, JHEP 0307 (2003) 
001. 



[7] J. Alwall et al, Eur. Phys. J. C 53, 473 (2008). 
[8] J. Alwall et al, JHEP 0709 (2007) 028. 
[9] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605 
(2006) 026. 

[10] E. Richter- Was,|arXiv:hep-ph/0207355| 
[11] M. Cacciari, arXiv:hep-ph/0607071 
[12] M. Cacciari, G. P. Salam and G. Soyez, JHEP 0804 
(2008) 063 

[13] M. M. Nojiri, Y. Shimizu, S. Okada and K. Kawagoe, 

JHEP 0806 (2008) 035. 
[14] M. M. Nojiri, K. Sakurai, Y. Shimizu and M. Takeuchi, 

JHEP 0810 (2008) 100. 
[15] J. Alwall, K. Hiramatsu, M. M. Nojiri, Y. Shimizu, work 

in progress 



